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We treat the problem of self-consistently interacting bosons in the presence of a finite (but 
macroscopic) potential well within a quasi-classical approximation for the normal component and 
Q ■ the order parameter. We solve the equilibrium problem and show, that actually condensation 

t^j- ■ occurs in two steps. One already at low densities with Bose condensation only in the well 

and another one corresponding to the usual condensation in bulk. The peak and width of the 
distribution of trapped particles in the well display a distinct signature of the local condensation. 
QO ' A possible connection to recent experiments with excitons is discussed. 
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£j ' 1 introduction 



co : 

Since Bose-Einstein condensation (BEC) of trapped atoms was experimentally put in evidence [U [2] 

c$ • much theoretical work was devoted to the discussion of the various aspects of BEC in an external 

£3 ■ potential (see Refs. 0IHE1 for recent reviews). 

! In this paper we develop a systematic quasi-classical approach for the description of weakly 

Ch ' interacting bosons with condensate and non-condensate in equilibrium and non-equilibrium. In order 
O 

o 



X 



to describe Bose-Einstein condensation of particles in a potential well which contains many closely 
spaced energy levels, a quasi-classical approach is appropriate, because the detailed knowledge of 
the population of higher lying energy levels is neither of much interest and may cause an enormous 
numerical effort for its determination. Our attention is devoted mainly to finite potential wells, for 
which the trapping in the well is controlled by the bulk population. In a previous publication [6j 
the evolution of a non-interacting Bose gas coupled to a thermal reservoir towards equilibrium in 
the presence of a finite potential well with many energy levels was discussed. It was shown that 
above a critical density the lowest state will be macroscopically populated, i.e., the local density 
will go to infinity. In this paper we take into account the repulsion between the bosons which 
eliminates this unphysical behavior. We develop this theory within a quasi-classical approach and 
a self-consistent treatment of the repulsion. We formulate the equations describing the evolution 
towards equilibrium for the normal component of the non-condensed particles, as well as for the 
condensate order parameter. We discuss in detail the exact equilibrium solution of these equations. 
The resulting scenario is a Bose condensation in two stages: At a first (rather low) critical overall 
density a locally confined condensate appears in the well, its distribution broadens with increasing 
overall particle density and at a second (higher ) critical density one gets the usual bulk condensate. 

* On leave of absence from Institute of Physics and Nuclear Engineering, Bucharest 
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The spatial distribution of the trapped particles changes abruptly as the local condensation in the 
well occurs. The width of the distribution displays as a function of the overall particle density a 
minimum in the form of a cusp, while the slope of the peak changes discontinuously. 

The described scenario may have a connection to recent experiments on excitons in a CU2O 
semiconductor under mechanical stress, where a strong narrow luminescence from the center of the 
well was observed. 



2 Quasi-classical description of a Bose system 

The Wigner function depending on the momentum p and the coordinate x is defined as 

dye^(ip(x + -y,t) + ip{x- ^y,t)) , (1) 

where ip{x, t) is the second quantized wave function and (. . .) means averaging over a given ensemble. 
This implies the normalization 

All the averages of operators which are a sum of two operators 0(x, |V) = A(x) + B(-V) 
depending on coordinate and momentum, respectively, may be expressed as integrals {0(x, fV)) = 
J dx J (2^3 f(p, x, t)0(x,p). In the quasi-classical limit f(p,x) is real and positive. For the case of 
an interaction with a thermal bath and slowly varying potentials [HI El the following bosonic kinetic 
equation holds 

(3) 

where the positive transition rates W(p,p') satisfy the detailed balance relation 



W(p,p') = W{p',p)e^ {p2 - p ' 2) , (4) 

and U(x,t) is a given potential. Generally speaking one has an external potential U (x), which we 
choose to be time-independent (the given potential well) and an internal one, due to the interaction 
between the bosons through a potential v(x) (where x is the relative coordinate of the two particles) 
related to the particle density n(x) (Hartree self-energy) 

U(x, t) = Uq(x) + J dx'v (x — x')n(x') = U (x) + J dx'v(x-x') J ^ /(p, x' , i) . (5) 

The kinetic equation Eq.® has a single stable solution, corresponding to equilibrium 

feM = -/„2 1 \ > (6) 

where \i is the chemical potential and U eq (x) is determined by the self-consistency equation 
U eq (x) = U (x) + f dx'v(x - x ')n(x') = U (x) + f dx'v(x-x') f fi^js feqfa % ') ■ ( 7 ) 
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The above described equations are valid only in the absence of a condensation. As it is easy 
to understand, above a critical overall density ^p- the Bose distribution in 3D cannot provide a 
description of the whole particle density, and this is in conflict with the conservation of the particle 
number. It can be shown rigorously [10J that for super-critical conditions, where a condensation 
occurs, Eq.Q has no solutions with ordinary functions but only with distributions. 

An appropriate approach to avoid this problem is to consider a two-component treatment where 
f(p,x,t) describes only the non-condensate and a condensate is taken into account explicitly. 

The existence of such a phase transition with spontaneous symmetry breaking in the presence of 
a potential is of course an idealization, which admits that the distance between the energy levels is 
negligible. This admition is implicitely made by the quasiclassical approximation. A discussion of 
the meaning of BEC in systems with discrete spectrum is given in Ref . [TT] . 

In the presence of a condensate the s.c. potential is then given by the total density (non- 
condensate and condensate) 

n{ ^ t) = J (2^/® f ^) + IW^))| 2 , (8) 

where (ip(xt)) is the condensate amplitude (average of the second quantized wave function) and thus 
the two components are coupled. A further coupling will result from irreversible transitions from 
and to the condensate. 

In order to include the condensate one needs an equation for its amplitude. For the mathematical 
consistency of the whole model, however, the quasi-classical approximation for the condensate also 
has to be formulated. One may start from the Schrodinger equation with the self-consistent (s.c.) 
potential U (x, t) for the condensate ampitude 

zh^(x, t)) = (~V + U(x, t)) (il>& t)) . (9) 



dt ' \ 2m 

Neglecting the contribution of the non-condensate to the s.c. potential, this equation coincides with 
the Gross-Pitaevskii equation, but we do not neglect this contribution. 

One has to recall here, that the quasi-classical limit means h — > and therefore in this limit the 
kinetic energy term disappears (but as usual the h from the time derivation survives) and we get 

iH^(x,t)) = U(x,t)(^(x,t)) . (10) 

Next, as already mentioned, we have to include the transitions between the normal component 
and condensate in the Eqs. (|3|l()|) . 

In the kinetic equation Eq.(|HJ) we have to add the scattering rate due to the transitions from the 
state p to the p = condensate 

-W{p,0)f{p,x,t)\(iP{x,t))\ 2 - (11) 

The initial state has to be populated, hence the factor f(p,x,t). The final state factor contains the 
term due to stimulated transitions oc | (ip(x, t)) | 2 , because the term describing spontaneous transition 
is lost in the thermodynamic limit (see also Ref. [Hi]). Similarly the scattering rate due to transitions 
from the condensate to a state p is given by 

W(0,p)(l + f(p,x,t))\(iP(x,t))\ 2 . (12) 
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Formally, these terms result by adding to the distribution function f(p',x,t) in the collision term of 
Eq.|2) a condensate term (2TTh) 3 5{p')\('^(x, t)}\ 2 . 

The gain and loss rates for the order parameter amplitude have to be calculated by taking the 
initial and final state factors for the condensate to be one 

Y<V>(i?,f)> / j^(w(p,0)f(p,x,t)-W(0,p)(l + f(P,^t))) . (13) 

Such collision terms have been obtained [E] in the frame of the equation of motion method for a 
system of bosons weakly interacting with a thermostat, decoupling higher order correlations and 
taking the long-time Markovian limit (see also jllj). A more involved approach, within the frame of 
quantum kinetic, taking into account time-dependent energy renormalizations due to boson-boson 
interaction was described in Ref. [12] • Neglecting these energy corrections and taking the long-time 
limit one recovers the previous results. 

The resulting system of coupled kinetic equations which describe the quasi-classical evolution of 
interacting Bosons within the self-consistent approximation as 

{jt + ^? Vs ~ VgU{ ^ t)v ") f 5 1] = ~ S WW 3 ( w ® f )fif ' f ' t} (1 + /07 '' f ' ^ ,} ) 

- (W(p, 0)m x, t) - W(0, p)(l + f(p, x, t))) I (V(5f, t))\ 2 , (14) 

together with 

ih^(x, t)) = U(x, t) t)) + ^ t)) J (W(p, 0)f(p, x, t) - W(0,p)(l + f(p, x, t)) 



(15) 



and of course with the self-consistency condition 



U(x,t)=U (x) + j dx'v(x-x')Q ^^f(p,x',t) + \(ij(x',t))\^j . (16) 
The equation for the condensate density |(^(z))| emerges as 

||^(f,t))| 2 = |^(x,t))| 2 | ^(^Oj/K^j-^^H/Kx.t))). (17) 
One sees that the equations conserve the total average number of bosons 



/ 



dx\{i;(x,t))\ 2 + JdxJ -J^f(p,x,t) = (N). (18) 

These equations describe an irreversible evolution. Above a certain overall particle density (at 
a given thermostat temperature) any arbitrarily small, but non-vanishing initial condensate in all 
points x will evolve to its finite stable equilibrium (non-homogeneous) value, while below this density 
it will disappear. 

Because the h — > limit for the Schrodinger equation is equivalent to the infinite mass limit, in 
this quasi-classical approximation the condensate does not propagate, but will be created or destroyed 
locally. 

A theory, in which the quasi-classical approximation is used only for the non-condensate, runs 
into mathematical inconsistencies. 
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The theory we described takes into account the effect of the particle repulsion for the spatial 
distribution only. Anomalous pair correlations giving rise to a modified quasi-particle (Bogolyubov-) 
spectrum were ignored. 

As for any quasi-classical approach we expect the validity of the above described theory for slowly 
varying potentials where the distance between the energy levels is very small. A detailed discussion of 
these conditions is given in the last section, where the theory is tentatively applied to interpretation 
of experimental data. 

For a local interaction v(x) = w5(x) the self-consistency equation becomes also local 

U(x,t) = U (x)+w(\(i;(x,t))\ 2 + J j^f(p,x,tfj , (19) 

which facilitates its treatment. In what follows we shall restrict the discussion to the equilibrium 
problem with such a local interaction. 



3 Bose condensation 

We discuss here the equilibrium solution of Eqs. f)14M16|) considering a local interaction with w > 
(repulsion) at a given temperature for various overall densities n = ^p- . Here the thermodynamic 
limit is implicitly understood. 

In the absence of an equilibrium condensate \(tp(x)) eq \ 2 = 0, the equilibrium solution for the 
normal component is the Bose distribution Eq.©. Above a certain overall density n (see Eq. (fT8|) 
divided by the volume V ) a local condensate appears. The normal component should be further 
distributed according to the Bose distribution, which is a stationary solution of the l.h.s. of the 
kinetic equation Eq.(JHJ) and of the part of the collision term which does not involve the condensate. 
The equilibrium condensate must be such as to cancel the additional collision term 

I t))\ 2 J (W(p, 0)f(p, x, t) - W(0,p)(l + f(p, x, *))) (20) 

of both equations Eq.JTJJ and Eq. (|T7l) . Inserting the Bose function, the equilibrium condensate 
density no(x) together with the detailed balance relation Eq. (JH), this term becomes 

no(x) J jMjs w (?> 0)f eq (p, x) (1 - e«cM*H0) . (21) 

As a result, the equilibrium condensate density Uq{x) may be different from zero only, where U eq (x) — 
fi = 0. Due to the positivity of f eq (p, x), the chemical potential always must lie below the minimum of 
the potential /i < minU eq . Therefore the condensate can exist only in the minimum of the equilibrium 
self-consistent potential. This again corresponds to the behavior of particles with an infinite mass. 

The order parameter itself cannot be constant in time, but according to Eq. (fT5|) oscillates as 
(ip{x, t)) eq = e~ l>Mt ^n (x) and it is actually determined only up to an arbitrary coordinate-dependent 
phase factor, which is determined from the initial condition. In general the order parameter in 
equilibrium always oscillates with the chemical potential because the particle number operator does 
not commute with the macro-canonical density matrix in the case of the spontaneous symmetry 
breaking (see [T2]). 

Without repulsion, the chemical potential can increase with the overall density up to the minimum 
of the potential Uo(x). Further increase of the overall density n is possible only, if the local condensate 
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density n Q (x) itself increases as the volume. This means that the assumed slow spatial variation which 
was required in the derivation of the quasi-classical description, breaks down. However, we shall see 
that the situation changes substantially by the existence of a repulsion between the bosons, allowing 
a quasi-classical description of BEG. 

Let us discuss, what happens in the presence of a repulsion when one increases the chemical po- 
tential starting form below the minimum of the external potential Uo(x). This is of course equivalent 
to increasing the overall density n. 

The self-consistency condition for equilibrium with a local interaction is 

U eq (x) = U (x)+w (n (x) + J ^^g /eg^x)^ . (22) 

This equation has been extensively used in the equilibrium theory of Bose-Einstein condensation of 
atoms trapped in an oscillator potential ( see Ref.jS] and references therein), where similar consider- 
ations apply, although the nature of the trapping kinetics is different. 

We choose an attractive external potential Uo(x) < which vanishes outside a finite domain 
around x = 0, and has a minimum in x = . One may expect that the potential U eq (x) also 
will have a point of minimum or a minimum surface and will be constant outside the range of Uq 
. Generally speaking minU eg > minU , since the repulsion opposes the attractive external potential . 

For a repulsive interaction the increasing density increases also the minimum of U eq (x) and even- 
tually flattens the bottom of the potential. The number of locally condensed particles will still remain 
finite and therefore does not contribute to the overall density. Only for an overall density n greater 
than the critical density 

We^v^L (23) 

an overall condensate appears. 

Because outside the range of Uq{x) both the local density and the self-consistent potential are 
constant, we have for n < n c the following equation for the determination of the chemical potential 

1 n. (24) 



(2^) 3 /(^+^) -1 

This equation has a solution only for /i < wn c . Above this value an overall condensate has to appear. 

For simplicity let us consider a spherically symmetric situation (Uo(r) with Uo(r) = for r > R 
and a negative minimum in r = ). We consider the case where /i is above the minimum of Uo, 
but below wn c . Let us look at the shape of U eq {r) coming from large r . Here one has a spatially 
constant self-consistent potential 

U eq (r) = wn; (r > R) . (25) 
For r < R one gets a radius-dependent density and the self-consistency equation 

f dv 1 

U eq (r) = U (r) + w / J! p2 ; (r < r < R) (26) 

for U eq (x) — \i > giving rise to a negative U eq (r) varying monotonously. At a certain r = r 
one may get U eq (r ) = \i and the integral achieves its maximal value n c . Due to the condition 
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U eq (x) — n > all the points r < r must belong to the minimum of a monotonous U eq (r) and a 
condensate n (r) must emerge 

U eq {r) = U {r) + w(n c + n (r)); (r < r ) (27) 
in order to ensure that U eq (r) remains constant for r < r , i.e. , 

t4j(r) = C/ e9 (r + 0) = /i (r < r ) , (28) 

implying 

n (r) = - (C/ (r ) - U (r)) (r < r ) . (29) 
w 

When the overall density n reaches the critical value n c , i.e. p = wn c , the self-consistent potential 
is completely flat, r approaches R, while the local condensate density approaches —-U (r) . 

For overall densities above the critical one, one gets an overall condensate, while the local con- 
densate density is n — n c — ^Uo(r). 

This scenario is complete and self-consistent! 

However, local condensation shows non-analyticities only of local entities, but does not affect the 
analyticity of the total thermodynamic entities. Therefore according to the usual terminology it does 
not define a phase transition. 

We expect the following scenario in a self-consistent but quantum mechanical picture: In the 
presence of a repulsive interaction between the particles the chemical potential approaches the lowest 
s.c. level and pushes it up without actually touching it until the bound states disappear completely 
and the true bulk phase transition occurs. In this way a finite local density appears first, which can 
be higher than the critical density of the bulk. 

4 Numerical illustration 

For the following we choose a spherically symmetric, attractive (v < ) external potential 

tW=*;(l-(£) 2 ) • (30) 

For further discussions it is convenient to measure densities in units of the critical density n c , the 
energies in units of ksT (U = f3U eq , Uq = (3Uq , M. = f3jj, , V = f3v, the radius in units of R (p = r/R 
), and use the dimensionless interaction constant W = f3wn c . 

Then 

U (p)=V(l-p 2 ) (V<0) (31) 
and one has to solve the equation 

1 f°° fs 

U( P)= U (P) + W—^ ds et+ J_ M _ i (32) 

JO e s — 1 

in every point po < p < 1 on order to reconstruct the condensate density as describe above. Obviously 
the radius Po = ^ for which the self-consistent potential touches the chemical potential is defined by 
the solution of the equation 

M = V (1 - p%) + W (33) 
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which exists in the interval < po < 1 only for A4 > V + W . 

Since the transcendental equation Eq. ([32|) is local, for a numerical solution it is convenient to 
solve it in favor of Wo at a given U, thus performing a simple integration. The association to a certain 
radius is given by the explicit definition of Uq(p). 

We illustrate the above described scenario with some numerical examples. We choose V = — 5 
and W = 0.05 and discuss later, what happens for stronger coupling interaction potentials. 

In FigUl the potentials Uo(p) (dotted line) and U(p) (full line) are shown for a chemical potential 
M. = —5.2 just below the minimum V of the external potential, i.e., slightly below the onset of the 
local condensation. One sees that the self-consistent potential U(p) is pushed up. For the chemical 
potential of M. = —4.8 (corresponding to a super-critical density) just above the minimum V of 
the external potential, Figj2] shows that the self-consistent potential obtains a flat bottom at the 
chemical potential. This is a typical signature for the appearance of a condensate here in the range 
< p < 0.2. FigJSlshows the resulting peak in the density. It is remarkable, that the overall densities 
for the two cases n = 0.002in c and n = 0.0031n c do not differ very much, both being much lower, 
than the critical one, while the condensate density in the second case exceeds the critical density by 
a factor 4 . 

A further increase of the chemical potential M. (or equivalently of the overall density n) will lift 
the bottom of the self-consistent potential further as it may be seen on FigQJ while the corresponding 
condensate density in the minimum of the self-consistent potential reaches very high values and the 
spatial density distribution (shown in FigEJ) widens out. Actually it is already very close to — 
, which will be its value at n = n c . 

The general feature of the condensation scenario is seen most clearly in FigEH where we represented 
the density peak of the trapped particles and its half-width as functions of the overall density n (here 
with the choice W = 0.05 and V = —5). The height of the density peak increases with the overall 
density but its slope increases abruptly at the overall density of 0.0027n c . The half-width of the local 
distribution at the same overall density has a minimum resembling to a cusp. 

For a stronger interaction, corresponding to W = 0.5 the role of the local condensate is not so 
striking. It can enhance the density in the well maximally by a factor of 10, resulting in a less 
pronounced peak. Nevertheless, the density inside the well may reach very high values compared to 
the overall density. This is illustrated in Fig. 0for two positions of the chemical potential M. = —4.7 
and M. = —4.3 corresponding to situations without and with condensate. In both cases the overall 
density is much below the critical one (n = 0.0035n c and 0.0052n c ). 

5 Excitons in a local potential 

In recent experiments Naka and Nagasawa have observed after uniform illumination a strong 
luminescence in the center of a strain induced potential well in a CU2O sample at T = 2°K. This 
luminescence may be due to the condensation of para-excitons (having a Bohr radius a ex = 0.7 x 
10~ 7 cm ). The luminescence peak was concentrated in a small vicinity around the potential well 
minimum. The strain induced potential in their experiment is asymmetric, but may be approximated 
qualitatively by a parabolic potential well like the one we choose having a depth of v ~ —lmeV and 
a radius of r ~ 10~ 2 cm. This would correspond in our dimensionless parameters to V = —5. In 
an ideal oscillator potential with the same curvature one gets an energy spacing of the quantum 
mechanical levels of ftw ~ 10~ 4 meV and a very small radial extension of the ground state oscillator 
wave function of 5 x 10~ 5 cm for the Cu 2 exciton mass of m = 2.7m . Since v/ftw ^> 1 many energy 
levels are contained in the potential well and one may expect the validity of the quasi-classical 
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description. At this temperature one has in Cu 2 a critical density of n c ~ 10 17 cm~ 3 . The density 
of the excitons created by laser illumination has not been measured directly. 

On the other hand, the interaction constant w is not directly experimentally accessible. Tenta- 
tively one may consider the well-known potential between hydrogen atoms obtained by the Heitler- 
London theory of the hydrogen molecule [TJ] and replace the Rydberg energy and the electron mass, 
with the corresponding entities of the exciton (exciton Rydberg Er and exciton Bohr radius a# 

v ( r ) = ~~F7~v? I d ^ / d y* ( (0(2/1)0(2/2) ± 4>(yi + r)4>(y2 - r)) 



S(r] 



1 + - - jj-— ) mmm) ) (34) 



\yi -m + f] \f\ \yi + r\ \y 2 - f| 
where the normalization factor is 



S(r) 2 = l±Jdy 1 J dyz^mWm + ?)<p(y 2 )<p(y 2 - f) (35) 
The Is wave-functions <f>(x) are 

1 s| 

0(f) = 3 . (36) 

Actually Sugiura [H] evaluated analytically all the involved integrals. The positive sign corre- 
sponds to the symmetrical state, while the negative one corresponds to the anti-symmetrical state in 
the coordinates. There are also more sophisticated variational calculations for the bound molecular 
state (symmetrical state), however not for the repulsive (anti-symmetrical) case, where a variational 
approach is not justified. 

The para-excitons (singlet spin state) in Cu 2 are the lowest exciton state. Then antisymmetric 
wave functions in the coordinates for the identical particles should be constructed and the resulting 
force between the excitons will be repulsive as shown in FiglHl The Heitler-London potential decreases 

2r 

with the distance between the excitons exponentially as e °« and the local approximation should 
hold. The calculated value of the interaction constant w is 49. 26E R a B respectively. For Cu 2 this 
value leads to W ~ 1. This means that in Cu 2 the strong coupling case is realized. 

Because we had to identify excitons with point like bosons, it is questionable to consider a 
Heitler-London potential between the excitons at distances smaller than 2ob- Actually one should 
consider such quantitative estimates with some reservations. Therefore, one cannot make a reliable 
quantitative prediction for the trapping wells in Cu 2 0. However, the photoluminescence of the 
trapped excitons should display the characteristics of FigEl if a local Bose condensation occurs. 

Finally we want to draw the attention to similar observations [H| with excitons confined to 
an indirect GaAs/AlGaAl quantum well layer. Here also a strong localized luminescence has been 
observed due to trapping, but the system is rather two-dimensional. 

Actually all the considerations relating to Bose-Einstein condensation of excitons are subject to 
the condition, that the electron-hole pairs are mainly bound into excitons, which under experimental 
conditions may be too restrictive. 

In conclusion, we have formulated a quasi-classical condensation kinetics of interacting bosons 
which are in contact with a thermostat and trapped in a finite macroscopic potential well. The 
theory is based on the coupled equations Eqs IT4lfT6l for the non-condensate and condensate. The 
equilibrium solution has been discussed in detail. Due to the repulsive interaction, the trapping 
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potential becomes renormalized after the condensation has occurred, allowing the chemical potential 
to increase over the minimum of the well for an increasing overall particle density. In particular it 
has been shown that the spatial distribution of the trapped particles shows distinct changes as a 
local Bose Einstein condensation occurs. The possible relation of these results to recent experiments 
with trapped excitons has been discussed. 

One of the authors (A.M.B) thanks the Deutsche Forschungsgemeinschaft for the generous support 
allowing his stay at the University of Frankfurt. We appreciate private communications on their 
experiments by N. Nagasawa and D.S. Chemla. 
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Figure 1: External potential Uo(p) (dashed line) and self-consistent potential U(p) (full line) for a 
sub-critical condition M. = —5.2 {n = 0.0021n c ) with a potential depth V = —5, and an interaction 
constant W = 0.05. 
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Figure 2: External potential U (p) (dashed line) and self-consistent potential U(p) (full line) for 
a super-critical condition M. = —4.8 ( n = 0.0031n c ) with a potential depth V = —5, and an 
interaction constant W = 0.05. 
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Figure 3: The local density n(p)/n c for the sub- and super-critical overall densities n = 0.0021n c 
( dashed line) and n = 0.0031n c (full line), respectively with a potential depth V = —5, and an 
interaction constant W = 0.05 . 




D 

Figure 4: The external potential U (p) (dashed line) and the renormalized potential U(p) (full line) 
for the normalized chemical potential M. = —0.1, and the corresponding density ( n = 0.584n c ) 
with a potential depth V = —5, and an interaction constant W = 0.05. 
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Figure 5: The local density n{p)/n c for the chemical potential M. = —0.1 (corresponding density 
n = 0.584n c ) with a potential depth V = —5, and an interaction constant W = 0.05. 
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Figure 6: The dependence of the height (thick line) and halfwidth (thin line) of the density dis- 
tribution of the trapped particles on the overall density in the vicinity of the condensation point. 
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Figure 7: The local density n(p)/n c for the sub- and super-critical overall densities at n — 0.0035n c 
(dashed line) and n = 0.0052n c (full line) for the strong interaction case with a potential depth 
V = — 5 , and an interaction constant W = 0.5. 
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Figure 8: The Heitler-London potential for the symmetrical and anti-symmetrical configurations 
(lower respectively upper curve). 
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